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Abstract 



Context. The synthesis of water is one necessary step in the origin and development of life. It is 
believed that pristine water is formed and grows on the surface of icy dust grains in dark interstel- 
lar clouds. Until now, there has been no experimental evidence whether this scenario is feasible 
or not on an astrophysically relevant template and by hydrogen and oxygen atom reactions. 
Aims. We present here the first experimental evidence of water synthesis by such a process on a 
realistic grain surface analogue in dense clouds, i.e., amorphous water ice. 
Methods. Atomic beams of oxygen and deuterium are aimed at a porous water ice substrate 
(H2O) held at 10 K. Products are analyzed by the temperature-programmed desorption technique 
Results. We observe production of HDO and D2O, indicating that water is formed under condi- 
tions of the dense interstellar medium from hydrogen and oxygen atoms. This experiment opens 
up the field of a little explored complex chemistry that could occur on dust grains, believed to 
be the site where key processes lead to the molecular diversity and complexity observed in the 
Universe. 
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1. Introduction 



Water, the spring of life (IBrack 2002r) . is the most abundant molecule in biological systems, and it 



is almost certainly of extraterrestrial origin. Water has been detected, in gaseous or solid form, in 
numerous astrophysical environments such as planets, co mets, interstellar clouds and star formin g 



regions where strong maser emission can be also observed jEhrenfreund et al. 2003 



Dartois 2005). 



Amorphous water ice was directly detected in dark interstellar clouds through infra-red absorption 
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dLeger et al. 1979h . During the formation of stars deep inside molecular clouds, gas and dust be- 
come part of the infalling material feeding the central object. Part of this gas and dust grains, 
covered with icy mantles (mainly composed of water), ends up in the rotating disks surrounding 
young stars and forms the basic material from w hich icy planetesima ls and later planets, together 
with comets in the external regions, are formed dvan Dishoeck 2004 ). W hile the means of deliv- 



ery of water to Earth remain a subject of debate d Morbidelli et al. 20001) . the synthesis of water 
in the Universe is a fundamental link in establishing our origins. Water molecule formation in 
the gas phase is not effi cient enough to reproduce the observed abundances in dark clouds, espe- 



cially in its solid form dParise et al. 2005 



Ceccarelli et al. 20071) . Therefore water ice must form 



directly on the cold interstellar grains and not as a condensate after being formed in the gas phase. 
A complet e review of th e processes involved both in the gas and solid phase has been recently 
publis hed ( Tielens 2005 ). It was suggested many yea rs ago that interstellar dust grains act as cat- 
alysts dOort & van de Hulst 1946c Ivan de Hulst 1946r) . Starting from simple atoms or molecules 
such as H, O, C, N, CO, grains are believed to be chemical nanofactories on which more com- 
plex molecules are synthesized leading eventually to prebiotic species produced concurrently by 
surface reactions and by UV photons an d cosmic rays irradiation, as already shown long ago 
dHagen et al. 1979t iPirronello et al. 1982 ). The most volatile species may be released in the gas 
phase upon formation dOarrod et al. 2007h . while the refractory ones remain on the grain surface, 
building up a so called "dirty icy mantle", a nd at least partially may be sputtered by the heavy com- 
ponent of cosmic rays d Johnson et al. 199 lb . Such mantles, having a typical thickness of a hundred 
monolayers, are mainly composed of water, the most abundant solid phase species in the Universe. 
Under dark cloud con ditions, except f or the very first monolayer that has to grow on bare silicate or 
carbonaceous grains dPapoular 20051) . most water molecules should be subsequently synthesized 
on a surface mainly composed of water. 

Chemical models in cluding water formation on grain surfaces were proposed years ago by 



dTielens & Hagen 19821) . They suggested that H2O formation would be initiated by H-atoms re- 



acting with O, O2 and O3, although the O3-1-H pathway was consi dered the most effective a nd O2 
would play more a catalytic role. Recent Monte Carlo simulations dCuppen and Herbst 20071) show 
that while the main route to water formation on cosmic dust grains in diffuse and translucent clouds 
is the reaction H + OH, in dense clouds the principal source of H2O is the reaction between H2 
and OH. This study also emphasizes the non-negligible contribution from the H + H2O2 reaction 
(H2O2 being a product of the H -H O2 pathway) and the unusual high abundance of reactants such 
as OH and O3 . Interestingly, another code by (author?) dParise 20041) proposed a water formation 
scheme where O3 molecules react with H- or D-atoms to form OH or OD, and subsequently the 
reaction H2 + OH/OD leads to H2O/HDO. It should be noted that this scheme was in part also 
constrained by the observed abund ances of deuterated s pecies. 

In previous laboratory works, dHiraoka et al. 19981) succeeded in producing water molecules 
from the reaction o f H- and O-atoms initially trapped in a N2O matrix. Very recently 
dMivauchi et al. 20081) investigated the reaction between cold H-atoms and an O2 ice at 10 K 
and demonstra ted the production of H2O2 and H2O molecules and estimated the efficiency of 
the reactions, dioppolo et al. 2008h did a similar experiment but with varying O2 substrate tem- 
peratures. They confirmed the production of H2O2 and H2O, made an estimate of the reac- 
tions efficiency and also drew conclusions upon the temperature dependence of the amount of 
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species produced. These two experiments dealt with the H2O production pathway in which O2 
is the species consumed to produce w ater as shown in preliminary experiments by our group 



dMomeni et al. 2007 



Matar et al. 2008) 



In the present study, the formation of water is studied for the first time using hydrogen and 
oxygen atoms interacting on the surface of an amorphous solid water (ASW) ice film, hence under 
conditions that are much more relevant to the interstellar medium. The aim of this first attempt to 
synthesize water under conditions close to those encountered in dense clouds was to investigate 
how water formation continues on the icy surfaces of cosmic grains and to give an estimate of the 
efficiency of the chemical path(s) actually active. 



2. Experimental procedures 



Experiments were performed with the FORMOLISM set-up dAmiaud et al. 2006h . In brief, a cop- 
per sample surface whose temperature can be controlled by computer in the 8-800 K range is 
maintained under UHV conditions; on it an amorphous solid water ice substrate, on which wa- 
ter formation is studied, is prepared in two step s. First, a 100-layer film of non-porous ASW ice 
is first deposited at 120 K ( IKimmel et al. 200 ih . then an overlayer of 10 layers of porous ASW 
ice is gro wn at 10 K. The unde rlying non-porous water film isolates the ice layer from the copper 
substrate (lEngquist et al. 1995b . This double ASW ice film is annealed to 90 K prior to each experi- 
mental run in order to avoid any further collapse of the pores between 10 and 80 K in the subsequent 
temperature-programmed desorption (TPD) experiments (see below). In fact, the species D2, O, O2 
used in this experiment will be thoroughly evaporated by 90 K and thus before any rearrangement 
in the porous structure of the ice template. 



The annealed water ice substrate is still porous dKimmel et al. 200 ll) and mimics an amorphous 



ice processed by UV and cosmic rays dPalumbo 2006h that is thought to c onstitute the icy mantl e 
of interstellar grains. By using an architecture with two separate channels (IPirronello et al. 1997b . 
each consisting of a triply differentially pumped beam line, two atomic beams of O and D are 
aimed at the ASW ice substrate held at 10 K. Atoms are produced by dissociation of O2 and D2 in 
microwave discharges. The dissociation efficiency of the oxygen beam is typically 40%, meaning 
that for 100 O2 molecules that initially feed the discharge, 60 O2 molecules and 80 O atoms will 
finally reach the cold target. Trace quantities of residual gases (i.e., CO, N2, CO2 and H2O) are 
present in the beam, although neither O3 nor deuterated compounds were detected. The D beam 
has a 60% dissociation efficiency, and no UV photons from the D2-discharge plasma can reach 
the water ice sample. The purity of the beams was checked with a rotatable quadrupole mass 
spectrometer (QMS) that intercepted the incoming beam. The O- and D-beams have a flux of 10'^ 
atoms/cm^/s and 5 x 10^^ atoms/cm^/s, respectively. After concurrent or sequential injection of the 
beams, a TPD is performed at a heating rate of 20 K/minute from 10 K to 200 K and simultaneously 
the desorbing species are monitored with the QMS. The signals of mass 19 (HDO) and mass 20 
(D2O) presented here are coiTected by subtracting the contribution of the H2O water substrate, 
which naturally includes some isotopes. 



Dulieu et al.: Water formation on interstellar dust grains 



5 



CD 
C 

'</) 
CO 

O 



2500 
2000 
1500 
1000 
500 



o HDO 

— D,0 

^— HDO in D only experiment 




80 100 120 140 160 

Surface temperature (K) 



180 



Figure 1. TPD profiles of Dj^O (thin line) and HD'^O (circles) after irradiation of the water sub- 
strate (H^^'O) with D atoms, '^O atoms and '^©2 molecules. Thick line: TPD profile of HD'^O, 
after 30 minutes of D irradiation of the water ice substrate. The water ice substrate is held at 10 K 
during all exposures. Traces are vertically shifted for clarity. 

3. Results and Discussion 

Several experiments have been performed. In the first one D-atoms were sent onto the water sur- 
face to confirm that hydrogen atoms do not react with water molecules in the substrate to pro- 
duce deuterated wate r molecules (see the solid thick line in Fig. [T]), as already mentioned by 



dNagaoka et al. 2005h . In addition, we determined that D2 molecules do not react with O-atoms 
nor with O2 molecules residing on the ASW ice surface. This, that may seem of secondary im- 
portance, is on the contrary a decisive result that proves that the water formation process requires 
hydrogen in atomic form. 

Fig. [T] shows the mass spectrum recorded during the TPD performed after the irradiation of 
the water substrate with D-atoms and '^O-atoms (and '^©2 molecules). Water formation clearly 
occurs as is testified by the presence of Dj^O and HD^^O peaks. These findings are not surprising 
although one could expect that D2O should be mainly formed, because D is used as precursor. 
But it is also know n that during the heating, isotope exchange between water molecules will occur 



jSmith et al. 19971) . For a very thin layer of deuterated water, a complete isotope exchange with the 
underlying H2O substrate molecules is expected around 150 K. In such a context, even if D2O is 
formed at 10 K, most of the deuterated water will desorb as HD^^O, as is observed in our experi- 
ments. 

Due to the mass spectrometric detection method of water molecules in the gas phase one could 
wonder whether molecules synthesized by the reaction of D-atoms and O-atoms and O2 were 
formed during the warm up of the whole ice substrate, a fact that would render the result not 
relevant to interstellar space. However, this is certainly not the case because reactions involving 
D-atoms must proceed at temperature below 20 K. At temperatures higher than that value the resi- 
dence time of H-/D-atoms becomes exceedingly small for reactions to occur. 

In order to confirm unequivocally that water was formed starting from the reactants deposited 
from the gas phase and provide quantitative results, we repeated the experiment presented above 
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Figure 2. HD'^O (open triangles) and Di'^^O (solid line) desorption profiles after ten minutes of 
simultaneous irradiation of the amorphous H2O ice substrate with '**0- and D-atoms. 

with the oxygen beam line fed with '**02 molecules. To be sure that the concentration on the sub- 
strate of D-atoms on the surface was always much higher than O-atoms and O2 molecules, we 
opened the oxygen beam by time-slices of 20 s every 2 minutes while the D-beam was constantly 
running and iiTadiating the target held at 10 K. Thus, the relative co ncentration on the surface con- 



sists of > 10 D-atoms per O-nucleus (in atomic or molecular form). (IMatar et al. 2008h have shown 



that D reacts very efficiently with O2, and that O2 disappears from the surface in a one D for on e O2 



fashion. In such a scenario, due to the very low mobility of O-atoms dTielens & Hagen 1982h and 
O2 molecules at the iiTadiation temperature and in the temperature range in which H- and D-atoms 
remain adsorbed on the substrate, it is reasonable to assume that O and O2 surface concentration 
reflects the proportions in which they were produced in the beam. O-atoms have therefore a very 
low probability to form O2 and to encounter O2 molecules to form O3, especially because they 
react with D with a higher probability. 

Fig.|2]shows a TPD mass spectrum obtained after a 10-minute dose (30 time-windows of 20 s 
every 2 minutes) of '^O-atoms and simultaneous D-atom irradiation of the ASW ice substrate at 
10 K. The two peaks at mass 21 and 22 can only be assigned to HD'^O and Dj^O and this clearly 
shows again that formation of water molecules has occurred on the amorphous H2O ice substrate, 
namely, a realistic analogue of cosmic dust surfaces in dense clouds. 

In Table[T|are shown the relative abundances (in terms of uncorrected ion counts) of the species 
detected during the TPD, and, in addition, the total number of oxygen nuclei ("^O atoms and twice 
the number of '**02 molecules) that have been exposed to the substrate during the irradiation phase. 
We have not performed here any coiTection (except the coiTection from the H2O substrate). As a 
preliminary remark, the signal at mass 23 is expected to be 0, and it serves as a good indicator of 
the level of noise. To evaluate the conversion efficiency of O-atoms into water molecules in our 
experiment, we may evaluate the ratio a between the sum of mass 20 ('^OD and H'^O ), mass 21 
(HD"*0), mass 22 (D^^O) and the total amount of "*0 nuclei: 

M20 + M21 +M22 

a = — = 0.76 

total ISO 
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Table 1. Surface area under TPD curves of different masses after the '^O + D experiment. Results 
of a statistical model are given as a comparison. 





Total '*'0 


Mass 19 


Mass 20 


Mass 21 


Mass 22 


Mass 23 


Area (Cts.K) 


15000+1000 


30000+2000 


10000±500 


1250+150 


250+70 


20 


Model 


15000 


31802 


8076 


1434 


188 






This ratio represents an upper limit to the formation of water molecules from O assuming 
that the contribution of Dj^O to the signal at mass 20 is negligible, as we have actually seen above 
(Fig.[T]). a, however, is not 1 . The missing '^O nuclei are likely to have partially desorbed during the 
reaction or ended up in the form '^OH radicals or even as H2^02 (and its deuterated isotopologues). 

The formation of water m olecules through the D + O2 pathway has already been stud ied and 



confirmed experimentally by (IMivauchi et al. 2008t lloppolo et al. 2008 



Matar et al. 2008h . In the 



present experiment, we can deduce a lower Umit for the formation efficiency for the D+O pathway. 
Actually 60% of the available '^O on the surface is contained in molecules. We can then see that, 
even if we assume that the molecular pathway has a 100% efficiency, it is less than a and a lower 
limit for water formation via the D+O pathway is about 0.16 (0.76 - 0.60). Therefore, if we 
define the formation efficiency via atoms as the ratio between the fraction of water molecules that 
certainly formed via the D + O route (0.16) and the fraction of O-atoms available (0.4), we obtain 
0.4 as a lower limit for the formation efficiency. On the other hand we can also deduce an upper 
limit. If we assume that O and O2 react with equal probability with D, we can estimate that 3 1 % 
(0.76 X 0.4) of the water we have formed is the product of the D + O pathway and the efficiency of 
water formation via atoms is equal to a. Considering both boundary conditions, we conclude that 
a reasonable estimate for the efficiency of water formation via atoms is about 0.5. 

In order to know if we can go deeper in the understanding of the different pathways of water 
formation, we have made a very simple statistical model which is able to reproduce correctly the 
data (see Table [1]). The model considers only statistical equilibrium. For sake of clarity we give 
here a simple example. Suppose there is an initial population of 90 H2O and 10 OD , the final 
equilibrium obtained by randomly mixing the isotopes is OD = 0.5, OH - 9.5 , H2O = 81.225, 
HDO = 8.55, D2O = 0.225. We have also included oxygen scrambling and calculated the statistical 
weight of all the final compounds. Finally, we sum the final products by masses taking into account 
QMS fractionation and are able to compare with the experimental results. 

The question is the following: is it possible that the OH + H reaction has a activation barrier as 
previously suggested in gas phase. In this case, '^OD could be mixed with H2O at higher tempera- 
ture during the isotope exchange? We first assume that the OD + D reaction is actually negligible 
compared to the O2 + D reaction, which leads to formation of OD and D2O in equal amounts. 
Therefore, at the end of the D irradiation the ratio a = ^^OD /D^^O is equal to (40+30 =70)/30 
because '^O atoms will form '**OD and '^02 molecules will form an equal amount of '^OD and 
D2'^0 (via '^02D2) . We assume only one adjustable parameter, the ratio fe='^0/'^0 of atoms that 
will exchange during the TPD and statistically equihbrate their populations via isotope exchanges. 

We find a best fit (Table [1]) for the only free parmeter b = 15. This means that the overlayer 
of products ('^OD and Dj'^O) is mixed with 15 layers of the H2O substrate and this is reasonably 
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consistent with previous results. (ISmith et al. 20001) found an efficient mixing of about 50 ML of 
both isotopes, although they used a lower heating ramp and a different amorphous ice. If one 
believes in our simple statistical model, it is also possible to explain our data starting with an high 
population of OD, and therefore that an activation barrier for the OD + D reaction exists at 10 K. 
Indeed, the isotope exchange blurs the chemical pathways, and the data provided here give not 
enough constraints to make a firm conclusion about which pathway is more efficient. None the 
less, on average, we have estimated that water formation via the D + O pathway is highly efficient 
(about 50%) and certainly is responsible for producing a significant fraction of water molecules. 
More laboratory investigations should be performed and we actually plan to do that by combining 
TPD and infra-red spectroscopy. 

The relevance of our results to astrochemistry is clear and gives strong experimental support 
to simulations of the formation and time scale of growth of water ice mantles in dense clouds. 
For the first time it is demonstrated on experimental grounds that water molecules can form on 
an amorphous water ice substrate under interstellar conditions (i.e., through surface reactions be- 
tween atomic hydrogen and atomic and molecular oxygen on ASW ice), allowing the growth of 
icy mantles that are observed in dense clouds. Also, in our experiments, as occurs on interstellar 
grains, the formation of water ice via reactions between hydrogen and oxygen atoms suffers from 
the competition between the formation of molecular hydrogen and molecular oxygen and probably 
this competition might reflect in the depth profile inside the ice layer and in the time evolution of 
the mantle itself. Several issues still need to be investigated: the efficiency at various coverages 
and with an abundance ratio between hydrogen and oxygen atoms closer to the interstellar one; 
the interaction of the atoms that land from the gas phase and those belonging to the surface; other 
chemical pathways that contribute to the formation of water and so on. Besides all these investiga- 
tions, of course, still to prove is how the first monolayer of water molecules may form efficiently 
on an amorphous silicate and/or and amorphous carbon layer If it can not form efficiently enough, 
it will be compulsory to accept that the very first monolayer has to be built by water molecules 
formed by gas-phase reactions and then accreted on the bare refractory grain surface; a fact of ma- 
jor interest inside dense clouds near the threshold of observability of icy mantles. These and other 
issues will be addressed in forthcoming papers. 



4. Conclusions 

We presented the first laboratory attempt to reproduce the formation of water molecules on of 
realistic space analogue of grain surfaces in dense molecular clouds. By exposing O- and D-atoms 
to an amorphous water ice substrate held at 10 K deuterated water molecules were formed with a 
high efficiency (~ 50 %). 
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